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The iron-pnictide superconductors have a layered structureformed by stacks 
of FeAs planes from which the superconductivity originates. Given the multi- 
band and quasi three-dimensional [Ij (3D) electronic structure of these high- 
temperature superconductors, knowledge of the quasi-3D superconducting (SC) 
gap is essential for understanding the superconducting mechanism. By using the 
fc^-capability of angle-resolved photoemission, we completely determined the SC 
gap on all five Fermi surfaces (FSs) in three dimensions on Bao.6Ko.4Fe2As2 sam- 
ples. We found a marked fc^ dispersion of the SC gap, which can derive only 
from interlayer pairing. Remarkably, the SC energy gaps can be described by 
a single 3D gap function with two energy scales characterizing the strengths of 
intralayer Ai and interlayer A2 pairing. The anisotropy ratio A2/A1, determined 
from the gap function, is close to the c-axis anisotropy ratio of the magnetic ex- 
change coupling Jc/ Jab in the parent compound [2j. The ubiquitous gap function 
for all the 3D FSs reveals that pairing is short-ranged and strongly constrain 
the possible pairing force in the pnictides. A suitable candidate could arise from 
short-range antiferromagnetic fiuctuations. 

Angle-resolved photoemision spectroscopy (ARPES) has played an important role in re- 
vealing the electronic structure of the pnictides. These measurements have typically been 
carried out at a fixed incident photon energy {hv) and varying incident angles that map out 
the planar band dispersion as a function of kx and ky. Thus far, four Fermi-surface (FS) 
sheets have been observed with two hole pockets centred around the F (0, 0) point and two 
electron pockets around the M (tt, 0) point in the unfolded two-dimensional (2D) Brillouin 
zone. Below the superconducting transition temperature Tc, nodeless superconducting (SC) 
gaps open everywhere on the FS sheets [3]-[6], pointing to a pairing order parameter with 
an s-wave symmetry in the a-h plane, in agreement with a number of theoretical results 
[71-[TT]. However there are other experiments which have indicated possible nodes in the 
superconducting gap of some pnictides, either line nodes in the a-h plane or nodes along 
the c-axis [T2H14]. It is well known that upon tuning the incident photon energy hv^ the 
allowed direct transitions will shift in energy and consequently in the momentum perpen- 
dicular to the a-h plane (fc^) which enables the determination of the electronic dispersion 
along the c-axis. In the free electron final state approximation, the conversion is given by 
kz = ^/2m[{hu — (f) — Eb) cos^ 6 + Vo]//i, where Vq is an experimentally determined inner 
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potential [l5]. Several earlier ARPES studies have used this A;2-resolving capability to probe 
the 3D dispersion of the normal state electronic structure in BaFe2As2-based pnictides, and 
found a large kz band dispersion in the orthorhombic phase in the vicinity of the parent 
compound [T6H20]. The 3D nature of the superconducting gap, which is critically important 
to the understanding of pnictide superconductivity, has yet to be studied. 

Fig. la shows the spectral intensity measured at 10 K at a photon energy of 46 eV 
which corresponds to the reduced = (see discussion below), plotted as a function of the 
binding energy and in-plane momentum along the F-X direction. Two hole-like bands are 
clearly observed, corresponding to the a (inner) and /3 (outer) bands respectively, in both 
the spectral intensity plot and the second derivative intensity plot (Fig. lb). The energy 
distribution curves (EDCs) show that the a band disperses towards the Fermi level (Ep) 
and bends back (Fig. Ic) as a result of the opening of the superconducting gap of ^12 meV, 
consistent with earlier ARPES experiments 0, ^ . When the photon energy is tuned to 32 
eV (fc^ = tt), in addition to the a and /3 bands observed around the zone centre (the Z 
point), a third hole-like band (Fig. Id) emerges between the a and /3 bands, while the band 
calculations [211423] have predicted the existence of the third hole-like band in the pnictides, 
mostly for the 1111 pnictides, there is no accurate prediction of the observed kz dispersion 
of this band for the 122 compound. We label this new band as the band, whose full 
dispersion and Fermi surface will be described in a separate paper. 

Comparing Fig. lb and Fig. Id quantitatively, it is clear that the a band dispersion 
undergoes a parallel shifts in binding energy under different photon energies, as shown in 
Fig. 2a. Moreover, this modulation exhibits a well-defined periodicity in the photon energy, 
indicating that the excess energy is carried by the coherent interlayer quasiparticle (QP) 
tunneling with well defined momentum kz- Converting photon energy into momentum kz 
with an inner potential of 15 eV [16j, we find that the periodicity in kz is remarkably close to 
the one expected from the lattice spacing between the adjacent Fe layers, i.e. 27t/c^ = 0.951 
A~-^, with =c/2 = 6.6 A (due to bilayer FeAs in Bao.6Ko.4Fe2As2) [2l]. To determine the kz 
dispersion of the a band, we stay at a fixed in-plane momentum such that the binding energy 
varies in a region sufficiently away from the SC gap, and plot the QP dispersion (Fig. 2b) 
as a function of kz as shown in Fig. 2c. The dispersion can be described remarkably well by 

-^soi^x) ky^ kz) = E2j^{kx^ ky) + 2t±co8kz (1) 



3 



with an interlayer hopping amphtude t± 2.3 meV in the binding energy range around 
20 meV (along the F-X direction). Here and thereafter for notational simphcity, k^^ky^kz 
carries the units of 1/a, 1/6, respectively. The kz dispersion results in the warping of 
the a FS sheet along kz- One important implication of the observed kz dispersion is that 
this band cannot be a surface state which would have no kz dispersion. To describe the 
underlying FS quantitatively, we extrapolate the fitted dispersion given in Eq. (1) to the 
Fermi level and display the in-plane Fermi wave vector [kp) along F-X for the a, (5 and 
the electron-like (7/^) bands as a function of kz in Figs. 2d— f, respectively. The FS area 
variation (defined as 5 in (1 ± 5)Afs) in the a-h plane are ^10%, 4% and 1% for the a, (3 
and 7/(^ bands, respectively, with the same periodicity along A;^. 

Having established the A;^-dispersion of the quasi-3D electronic structure, we turn to the 
kz dependence of the SC gaps on diflFerent FS sheets obtained using many photon energies. 
Fig. 3a shows a collection of the EDCs at the Fermi crossings of the a band for different 
photon energies /iz/ = 30 — 60 eV. Appreciable gap variations are clearly visible in the 
symmetrised EDCs in Fig. 3b. The extracted SC gap values (defined as Ag^) at both left 
and right Fermi crossings are plotted in Fig. 3c as a function of the photon energy (left 
axis) and c-axis momentum kz (right axis). Remarkably, shows rather large periodic 
variations from ~9 meV to ^13 meV, then back to ^ 9 meV as kz moves from Z to F and back 
to Z. Similar kz dependence of the smaller SC gap on the (3 band is also observed, varying 
from ^5 meV to ^7 meV, as shown later in Fig. 4. However, the kz variation of the SC gap 
on the electron-like (7/^) FS sheets is much smaller, as indicated in Figs. 3d— f where the 
Bogoliubov QP peak is shown to situate at a nearly constant energy of ^11.5 meV when kz 
varies from F to Z. We note that an appreciable A;^-dependence of pairing strength associated 
with 3D band structure has been predicted in this material [25j, although the predicted 
certain in-plane anisotropy and gap nodes have not been observed in our experiments. 

It is interesting to note that the observed gap values (Figs. 3c and f) and the FS warping 
(Figs. 2d and f) along kz are anti-correlated, that is, when the in-plane FS area at a fixed kz 
becomes larger (smaller) the SC gap becomes smaller (larger) on the same plane. At the first 
glance, one might conclude that the gap variation along kz originates predominantly from 
the tunneling induced FS warping, that is, the kz dependence of the in-plane Fermi vector 
fci?, as is implied by the simplest form of an in-plane 5± gap function A^^ = Aq cos kx cos ky. 
However, this turns out to be not the case. The < 10% change in the in-plane Fermi vector is 
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too small to account for the large, near 40% gap variations because of the small "gap velocity" 
in the pnictides. This is clearly seen from Fig. 1 where the near doubling of the Fermi vector 
in going from the a-FS to the /3-FS only results in the gap change from 12 meV to 6 meV. 
To illustrate this point further, we plot in Fig. 4a the measured gap values as a function 
of I co8kxC08ky\ at the Fermi points for the a, a', /3, and the electron (j/S) bands. Notice 
that although the average gap value follows this 2D form, the marked deviations induced 
by the dispersion could indicate that pairing is not purely two-dimensional and there is 
an additional driving force, namely, the pairing between the layers that is predominantly 
responsible for the gap dispersion with k^. 

For an anisotropic layered superconductor with interlayer coupling, we adopt a simple 
form of the gap function, 

^3D{kx, ky, k^) = A2D{kx, ky){l + rjcosk^). (2) 

This is a direct generalization of the expression for BCS superconductors with an isotropic 
in-plane gap function [26j. Since the FS warping is rather small, A2D{kx^ ky) is approximately 
a constant. We expect this expression to be a good approximation, where 77 is a measure 
of the ratio of the interlayer to in-plane pairing strength. In Fig. 4b, we plot the measured 
A'^D on different FS sheets as a function of cosfc^. To increase the accuracy, we only used 
the lower photon energy part of the gap dispersion and averaged over the left and right 
crossings. Figure 4b shows that Eq. (2) fits the data rather well, with the anisotropy ratio 
Tj ^ 0.17, 0.13 and -0.01 for the a, /3 and '^/S bands, respectively (the gap values of the 
a' band are only resolvable near k^ = tt). It is interesting to note that the values of rj 
for the a and /3 bands are consistent with the ratio of the exchange coupling Jc/ Jab in the 
magnetically ordered parent compound extracted from the spin wave dispersion measured 
by neutron scattering [2j, where Jc ^ 5 meV is the interlayer coupling and Jab ^ 30 meV 
is the next nearest neighbour coupling in the Fe-plane. Note that fc-averaged Jc would 
have most contributions from the bands with appreciate kz dispersion, such as a, a' and /3 
hole-like bands. 

The observation of the cosine dependence of the SC gap on kx^ ky and k^ indicates that the 
gap function refiects the lattice symmetry and that the predominant pairing is short-ranged 
in real space. Under such an assumption and taking into account the lattice symmetry, 
the leading terms of a generalized 5-wave gap function can be written as Aicoskxcosky + 
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^{coskx + cosky)co8kz + As{co8kxC08ky)co8kz' It is possible that the gap parameters of 
Ai, A2 and A3 have some band/orbital dependence. For simphcity, we choose band/orbital 
independent parameters to fit the data. We found that A3 <C A2 as a result of the vanish- 
ingly small r] of the ^/S electron-like FS around (tt, 0). Remarkably, the remaining terms 
A 1 cos A; a; cos A;^ + ^(cosA;^; + co8ky)co8kz fit all the gap values on the different FS sheets 
quite well, with Ai = 12.3 meV, and A2 = 2.07 meV as shown in Fig. 4c. The ratio of 
A2/A1 is nearly the same as the ratio of Jc/Jab- We have also checked that the SC gap on 
each observed FS sheet along different in-plane directions fits well to this 3D gap function. 

Our finding of a single three-dimensional superconducting gap function for all five different 
FS pockets indicates that there are only two dominate pairing energy scales, one in-plane 
and one out-of-plane, in pnictide superconductors. It points to a common origin for the 
pairing strengths on all the observed FS sheets, independent of their different density of 
states. Moreover, because this gap function can be obtained by decoupling the 3D next 
nearest neighbour antiferromagnetic exchange couplings within the pairing channel, our 
results are consistent with short-ranged antiferromagnetic fiuctuations induced pairing in 
the iron-pnictide superconductor P [IHl Elj . 

Methods 

We have performed systematic photon energy dependent ARPES measurements in the 
superconducting state of the optimally hole-doped Bao.6Ko.4Fe2As2 superconductors (Tc = 
37 K). High quality single crystals used in our study were grown by the fiux method [28j. 
High-resolution ARPES measurements were conduced at the SIS beamline of the Swiss 
Light Source. The photon energy used in the experiments is between 20 eV and 110 eV with 
different circular and linear polarization. The energy resolution is 8 - 20 meV depending on 
the photon energy, and the momentum resolution is below 0.02 Samples were cleaved 
in situ and measured at 10 K in a working vacuum better than 5xl0~^^ Torr. The mirror- 
like sample surface was found to be stable without obvious degradation during a typical 
measurement period of 24 hours. Many samples have been measured and reproducible 
results have been obtained in these samples. 
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FIG. 1: Band dispersion of superconducting Bao.6Ko.4Fe2As2. a, ARPES spectral intensity 
measured at 10 K plotted in false color scale as a function of the in-plane momentum (fc||) and 
binding energy along F-X by using 46-eV photons which corresponds to kz = 0. Two hole-like 
bands {a (inner) and /? (outer)) are observed, b, Second derivative of the spectral intensity plot as 
shown in panel a. c, A set of EDCs within the E-k range indicated by the red rectangle in panel 
a. The red EDO is at A;^ of the a band, d, Second derivative plot of the dispersion along Z-R 
{kz — tt) measured by using 32-eV photons. Three hole-like bands {a (inner), a' (middle) and /? 
(outer)) are observed. 
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FIG. 2: kz dispersion of quasiparticles and Fermi surface warping, a, Comparison of the 
a band dispersions at A;^ = (red dots) and kz — tv (black dots) extracted from the peak position 
of momentum distribution curves along F-X. The solid lines are linear fits to the corresponding 
dispersions. The green dots denote the binding energies of the quasiparticle peak of the a band as 
shown in panel b measured at a fixed fcy (slightly below kp) by using different photon energies (or 
different kz). b, Corresponding EDCs to the green dots shown in panel a. The black vertical bars 
indicate the binding energies of the QP peaks. The red (blue) EDC corresponds to A;^ = {kz = 
tt). Note that the "shoulder" at the higher binding energy (^20 meV) of the EDC measured at hu 
= 58 eV is from the band, c, kz dispersion of the a band extracted from panels a and b. The 
solid curve is the fit using cosA;^. d, kz dependence of the in-plane k^ of the a-FS parallel to F-X. 
The solid curve is the cosA;^ fit. e and f, Same as panel d, but for the /3-FS and electron-hke (7/^) 
FS. 
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FIG. 3: fc^ dependence of the superconducting gaps, a, Photon-energy-dependent EDCs 
measured at the left on the a-FS along T-X or its parallel directions with different kz- Red 
(blue) EDCs correspond to A;^ = {kz = tt). The black blocks indicate the binding energies of 
the coherent peaks, b, Corresponding symmetrised EDCs of the ones shown in panel a. The 
dashed line at 12 meV is a guide to eyes for viewing the variation of the SC gap at different hu. c, 
Extracted values of the SC gap (defined as the half value of peak-to-peak positions in symmetrised 
EDCs) on the a-FS at different photon energies. The dots (squares) are obtained from left (right) 
side of A; on the a-FS. d-f, Same as panels a-c but on the electron-like FSs (7/^) along F-M or 
its parallel directions. The dashed line in e is at 11.5 meV to guide eyes for viewing the SC gap at 
different hu. 
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4c = ^2D (1+rjcoskJ 



— A2d" = 11-1 meV, 11^ = 0.172 

— A^d'^ =6.36meV,f]f^ = 0.131 

— A^^y^^ ^11.4meV, rf^ = -0.012 



A^^ = Afiosk^cosky +-^(cosk^+cosky)cosk^ 
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FIG. 4: 3D superconducting gap function A{kx^ ky, kz)- a, The SC gap values on 
the a-FS (red dots), /3-FS (blue dots), j/S-FS (green dots) and a'-FS (brown dots) as func- 
tions of I cos /c^; cos I . The black solid line is the gap function |A| — | cos /c^; cos /c^ | with 
Ao = 12.3 meV. b. Same as panel a but as functions of cosA;^. The solid lines are indepen- 
dent linear fits to the SC gaps on the different FS sheets using a generic 3D gap function 
\A{kx,ky,kz)\ — \A2D{kx,ky){l^r]coskz)\- c. Same as panel b but using a single 3D gap func- 
tion \A{kx^ky,kz)\ = I AicosA;a^cosA;y + ^(cosA;^^ + cosky)coskz\ with Ai = 12.3 meV and A2 = 
2.07 meV to fit ah the SC gaps. 
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